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SECONDARY PATTERN COMPUTATION OF AN OFFSET REFLECTOR ANTENNA 


Roberto J. Acosta 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

Reflector antennas are widely used In communications satellite systems 
because they provide high gain at low cost. In analyzing reflector antennas 
the computation of the secondary pattern Is the main concern. A computer pro- 
gram for calculating the secondary pattern of an offset reflector has been 
developed and Implemented at the NASA Lewis Research Center. The theoretical 
foundation for this program Is based on the use of geometrical optics to des- 
cribe the fields from the feed to the reflector surface and to the aperture 
plane. The resulting aperture field distribution Is then transformed to the 
far-fleld zone by the fast Fourier transform algorithm. Comparing this tech- 
nique with other well-known techniques (the geometrical theory of diffraction, 
physical optics ( Jacobl-Bessel) , etc.) shows good agreement for large (diameter 
of 100X or greater) reflector antennas. 

INTRODUCTION 

The accurate prediction of radiation characteristics for a microwave 
antenna Is essential In designing antenna systems. Antenna radiation charac- 
teristics such as beam width, gain, aperture efficiency, side-lobe level, and 
cross polarization are used In analyzing and designing advanced antenna sys- 
tems. The aperture field method (ref. 1) described In this report Is one of 
several methods for predicting antenna performance characteristics. The method 
assumes that the tangential electric field on a planar surface In front of the 
antenna reflector Is known. This aperture field distribution Is then trans- 
formed to the far-fleld zone of the reflector by using a two-dimensional 
Fourier transform. 

To compute the tangential electric fields on the aperture surface, the 
geometrical optics method (ref. 2) Is used. In this technique the energy com- 
ing from the feed and reradiating from the reflector surface Is characterized 
by an astigmatic tube of rays (ref. 3) that allows calculation of the amplitude 
and polarization of the electric field. This report also presents a method 
for computing the secondary pattern by using the aperture field method and a 
fast Fourier transform (ref. 4). Computation of the tangential electric field 
In the aperture plane Is briefly described. A reflector configuration (fig. 1) 
Is analyzed, and the results are compared with other well-known computational 
techniques. A description and a copy of the program are Included In appendixes 
A and B, respectively. 

The author wishes to thank Dr. Shung-Wu Lee and Dr. Peter T.C. Lam, of 
the University of Illinois, for their cooperation and recommendations. 





APERTURE FIELD METHOD APPROACH 
Description of Problem 

The geometry of the problem under consideration Is shown In figures 2 
and 3. A reflector Er Is Illuminated by the energy from a feed source at 
Fi . The problem Is defined as (1) to determine the tangential electric field at 
an observation point on the aperture grid F 2 , as shown In figure 2 , and ( 2 ) 
to determine the secondary pattern from the near-field distribution by using 
the fast Fourier transform, as shown In figure 3. 

Incident Electric Field on Reflector 


The radiated electric field from the feed antenna has the asymptotic form 
given by equation ( 1 ) : 


-Jkr . 

t - v ?<».»> I -■ ( 

where f* - Is the active element pattern, k = 2 ir/\ Is the wave number, 
and r Is the distance from the source to the reflection point. The vector 
function In equation ( 1 ) can be approximated by equation ( 2 ). 


? (0,9) £ eU E (©)(ae^ cos 9 + b sin 9) + $U H (e)(b cos 9 - ae^ sin 9) 

where Uf(0) Is the E-plane active pattern, Uh(0) Is the H-plane active 
pattern, and a, b, and 9 are polarization parameters: 

The various feed polarizations are described In the following table: 

a b 9 

Linear x 100 

Linear y 010 

Right-hand circular polarized (RHCP) 1/^2* 1/^T ir/2 

Left-hand circular polarized (LHCP) 1 /^ 2 " -«/2 

Typically these element patterns can be approximated by (cos 6)9; that Is, 


( 2 ) 


u E (e) = (cos 0 ) 


U H ( 0 ) = (cos 0 ) 


(3a) 

(3b) 


If equations (3a) and (3b) are used to represent the active element patterns, 
the power radiated (ref. 5) by this source Is given by equation (4). 
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Two coordinate systems (fig. 4) are used, the feed coordinate system 
(x F ,y f ,z F ) and the reflector system (x,y,z). These two coordinate 
systems are related by Eulerlan angles (tt,B,y) Illustrated In figure 5. 
Reference 6 provides a good description of these angles, which transform the 
Incident field on the feed coordinate system Into the reflector coordinate 
system. 


Reflected Electric Field 

For a given feed point F-| and an observation point F 2 (fig- 2), a reflec- 
tion point 0 r may exist on the reflector ip. This point Is called the 
specular point. This type of ray reflection satisfies Snell's law of reflec- 
tion. Reference 7 describes a search procedure to obtain the specular point. 

The reflected field at the aperture point F 2 Is given by equation (5). 


t[ F2) - OF • e Jkd2 [2(n • - t 1 ] j . (5) 


where d 2 Is the distance from 0 r to F 2 . 

Equation (5) Is given In terms of the Incident field at the reflection 

point 0 , the surface unit normal n of the reflector at 0 /’and a 
divergence factor OF. The divergence' factor In equation (5) Is given by 


( 6 ) 


where (r|, R^) are the principal radii of curvature of the reflected wavefront 

passing through 0 r . Reference 7 gives an expression for these parameters 
as a function of the principal radii of curvature of a parabolic reflector. 



Secondary Pattern 

. From the Hygens theorem (ref. 1) a solution for the far-fleld zone may be 
obtained If the tangential fields at (fig. 3) are known. The aperture 
plane Is taken to be perpendicular to the z-axls. 

The field In the far-fleld zone Is given In terms of the tangential 
electric field at Z A , which Is denoted by t (ref. 8). Equation (7) relates 

A Q 

the aperture field ? to the far-fleld zone electric field. 

3 
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Ee = Fy cos ip cos e - F x sin <p cos e 

(7a) 

E<p * F x cos «p + Fy sin <p 

(7b) 

^//V k,ux * vy > 

£ A 

(7c) 

u = k sin e cos <p 

(7d) 

v s k sin e sin <p 

(7e) 


and e and <p are the spherical coordinates of far-fleld point. The Integral 
In equation (7c) Is a double Fourier transform. The computer program developed 
uses a fast Fourier transform (FFT) algorithm to accomplish this task. Refer- 
ences 9 and 10 describe In detail the use of the FFT algorithm to solve for 
the double Fourier transform. « 

The aperture field theory used to determine the secondary pattern Is exact 
If the tangential fields are known everywhere on the aperture plane r A . When 
employing the FFT, 1/^ Is truncated. To minimize the amount of computer time, 
z a should be as small as possible while capturing almost all of the field. 

The number of grid points should satisfy the Nyqulst (sampling) theorem. The 
secondary pattern can be reconstructed from a discrete set of near-field values 
If the spacing Is X/2 or less. Reference 11 describes sample spacing for ^ 
high-gain reflector antennas. 

This method of calculating the secondary pattern^ls accurate In cases v ^ 
where the antenna diameter Is of the order of 100X or more. If the antenna \ 
diameter Is less than 100X, the accuracy Is reduced, specifically In the slde-y 
lobe region. This approach was originally used by fr uf -f man - and C n 4 we l T 7 reTTTJ 
and developed extensively by many others (e.g., Hwang, et al., ref. 12). The 
reflector configuration described In figure 1 was analyzed by using various 
methods (physical optics ( Jacobl-Bessel) (refs. 13 and 14), the geometrical 
theory of diffraction (ref. 15), and geometrical optics). The calculated pat- 
terns and directivities are shown In figures 6 and 7. Figure 6 shows an 
on-focus case; the directivity and the far-fleld pattern are In very good 
agreement. Figure 7 shows an off-focus (or scan) case; the agreement Is still 
very good. The computer program given In appendix B was used to analyze this 
configuration. ?lin tti irtlv-Uy . Appendix A 

contains a detailed description (users guide) of all the Input parameters. 
Including those for the configuration shown In figure 1. 


Directivity 

The far-fleld zone Is usually divided Into two orthogonal polarizations. 
Following Ludwig's definition 3 (ref. 16), the following unitary polarization 
vectors are Introduced: 

R = e(ae^ cos (p + b sin <p) + <p(-ae^ sin <p + b cos <p) (8a) 

C = e(ae - ^ sin «p - b cos <p) + $(-ae"^ cos <p + b sin <p) (8b) 
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If the secondary pattern can be expressed as 


-Jkr 

f - V t §t 9 * * E ,> <9) 

the reference-polarization expression for ? Is 

t • (R*)* (10a) 

and the cross-polarization expression Is 

t • (C*)* (10b) 

The directivity for the reference polarization Is defined by 

* 2 

4* IE • R| /Z n 

0 R (e,«p) = P * ( 11 a) 

rad 

Similarly directivity for the cross polarization Is defined by 

4* | t • C| 2 /Z Q 

D r (©,«p) p (lib) 

u p rad 

CONCLUDING REMARKS 

A method, called the aperture field technique, has been developed for 
calculating the secondary pattern of an offset reflector Illuminated by a feed 
with arbitrary polarization. By using the fast Fourier transform the far-zone 
electric field Is computed very efficiently. The results for the secondary 
pattern are In good agreement with those obtained by other well-known tech- 
niques. This method can be conveniently extended to secondary pattern compu- 
tation for multiple and numerically specified reflectors. 

The computer program based on the aperture field technique Is one of the 
main research tools used at the NASA Lewis Research Center for analyzing 
advanced antenna systems. 
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APPENDIX A 


DESCRIPTION OF PROGRAM 

A computer program was designed to calculate the antenna far-fleld 
performance characteristics. The method of analysis Is geometrical optics. 
The output (unit FT07) Is Interfaced with a plotting program to display the 
far-fleld values. Any of the available plotting routines In the computer 
library may be used. The main features In this program (SOURCE • RAP$11) 
are the calculation of directivity, feed efficiency, and plots of E or H 
far-fleld plane cuts. 


Input (Unit FT05) Parameters to RAP*11 


F 

FI 

F2 

DIA 

XF, YF, ZF 

DELTX 

DELTY 

FREQ 

QE 

QH 

POL 

ZETA 


focal length, m 

fraction of wavelength that defines surface sample spacing 
fraction of wavelength that defines aperture sample spacing 
diameter of reflector (If circular projection), m 
coordinates of center (origin) of feed coordinate system on reflector 
coordinate system 

offset of center of reflector relative to x-axIs, m 
offset of center of reflector relative to y-axls, m 
frequency 

E-plane cosine pattern exponent 
H-plane cosine pattern exponent 
1, Y-POL; 2, X-POL; 3, RHCP; 4, LHCP 

maximum far-fleld angle for which pattern will be calculated 


The Input Is described In figure 8. 


The Input for the case described In figures 1, 6, and 7 (to be read In 
FT05) Is as follows: 


F = 2.43 

FI = 0.5 

F2 = 0.5 

DIA = 2.754 

XF = YF = ZF = 0.0 

DELTX = -1.377 


DELTY = 0.434 
FREQ = 11.74E9 
QE = 3.2 
QH = 2.6 
POL = 3 
ZETA = 4.0 
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APPENDIX B 


COMPUTER PROGRAM 


0000100 

0000200 

0000300 

0000350 

0000400 

0000410 

0000420 

0000425 

0000500 

0000501 

0000502 

0000503 

0000504 

0000505 

0000506 

0000507 

0000508 

0000509 

0000510 

0000550 

0000600 


0000700 
0000701 
0000710 
0000801 
0000802 
0000900 
000100 0 
0001100 
0001200 
0001300 
0001400 
0001500 
0 00 16 00 
00017 00 
0001710 
0001800 
0001900 
0002000 
0002100 
0002110 
0002120 
0002130 
0002131 
0002140 
0002150 
0002151 
0002160 
0002161 
0002162 
0002165 
0002166 
0002168 
0002169 
0002170 
0002171 
0002175 
0002176 
0002180 


C*hk#h*THIS PROGRAM IS REFLECTOR ANALYSIS THIS THE MOST GENERAL VERSION 
CxxxxxxTHE METHOD OF ANALYSIS IS APERTURE INTEGRATION. 

DIMENSION PNUMC 50 0 * 500 ) 

DIMENSION X(500),Y(500),ERX(500,500),ERY(500,500) 

DIMENSION PH( 500,500 ) , EXRC5000 ) , EXI (5000) ,EYR(5000) ,EYI (5000) 
DIMENSION FIELDXC 10 00),FIELDYC 1000 ) ,FIELD( 1000) ,ANGC 1000) 
DIMENSION RTNARRC2) 

CxxxxxiNPUT PARAMETERS TO THIS PROGRAM***************************** 
CxxxxxF • FOCAL LENTGTH M 

C*xx**Fl : FRACTION OF WAVELENTGH (SURFACE SAMPLE SPACING) 

C*x*x*F2 : FRACTION OF WAVELENTGH (APERTURE SAMPLE SPACING) 

C*xxx*DIA 5 DIAMETER M 

C*xxx*XF,YF,ZF : ORIGIN OF FEED COORDINATE SYSTEM. 

C*XX*XDELTX,DELTY s OFFSET X-AXIS # Y-AXIS RESP. 

C*xxx*FREQ • FREQUENCY HZ 
C**X*XQE,QH,Q : COSINE PATTERN EXPONENT 
C*****POL : 1 Y-POL, 2 X-POL, 3 RHCP, 4 LHCP 
C*x*x*ZET A : MAX FAR-FI ELD ANGLE. 

INTEGER*2 N1 

NAMELIST/INPUT/F,F1,F2,DIA,XF,YF,ZF,DELTX,DELTY,FREQ,Q,ZETA,QH, 
QE, POL 

READC5, INPUT) 

C*xxx*ALL CONSTANT PARAMETER IN THE PROGRAM*************************** 
PINT=0. 

WRITE(9,997) 

997 FORMAT ( 5X, 1 THIS IS RFLECTOR ANALYSIS PROGRAM. BY R.J. ACOSTA*) 

PI = 4 . *ATAN ( 1 . ) 

AWAVE=3E8/FREQ 
SSP=F1*AWAVE 
NSP= ( DI A/SSP ) +1 
SAP=F2*AWAVE 
NAP=(DIA*1. 5/SAP )+l 
RMAX=0 .5*DIA 
XC=DELTX+RMAX 
YC=DELTY+RMAX 

ZC=-(XC**2+YC**2)/(4.*F)+F 
XSMAX=XC+RMAX 
YSMAX=YC+RMAX 
XAMAX=XC+RMAX*1.5 
YAMAX=YC+RMAX*1 .5 
ALIM=.5*SAP 

FACT = ( (AWAVE**2)*(120 .#PI> > 

FACT = 1/FACT 

CxxxxxpOWER RADIATED BY THE FEED SOURCE . xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PRAD=60.#(2.*Q+1) 

PRAD=1/PRAD 

C*****THIS CALCULATES THE SUBSTENTED ANGLE BY REFLECTOR REL TO FEED***** 
XU=XC 

YU=YC+RMAX 

ZU = - (XU**2+YU**2 )/ ( 4 . XF) +F 
XL = XC 

YL =YC-RMAX 

ZL=-(XL**2+YL**2)/(4 .*F)+F 

A117=(XU-XF)X(XL-XF)+(YU-YF)X(YL-YF)+(ZU-ZF)X(ZL-ZF) 

A1 15=(XU-XF) **2+ ( YU-YF) **2+ (ZU-ZF) **2 
A115=SQRT(A115) 

A116=(XL-XF)**2+(YL-YF)**2+(ZL-ZF)x*2 
All 6 =SQRT ( A1 1 6 ) 

SI=ARC0S(A117/(A115*A116>) 
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0002181 

0002185 

0002186 

0002187 

0002188 

0002189 

0002190 

0002195 

0002196 

0002200 

0002210 

0002220 

0002230 

0002240 

0002250 

0002260 

0002300 

0002309 

0002310 

0002311 

0002315 

0002320 

0002325 

0002326 

0002327 

0002328 

0002330 

0002350 

0002360 

0002370 

0002371 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004800 

0004900 

0005000 

0005100 

0005200 

0005300 

0005400 

0005500 

0005600 


CxxxxxcALCULATION OF POWER INTERCEPTEDxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx . 
EXP=2.*Q+1 

Pll=-COS< .5*SI)XXEXP 
P12=P11/EXP+1/EXP 
PINT=P12/60 
PEFFR=PINT/PRAD 

TAPPER=20.XALOG10(COS(SIX.5)XXQ) 

WRITEC 9, 998 )PEFFR, TAPPER 

998 FORMAT (/////, 5X, * EFFICIENCY OF FEED 1 , F10 . 3 ,5X, ’TAPPER OF 
.FEED =’,E15.5) 

DO 400 1=1, NAP 
DO 410 J=1 , NAP 
PNUMC I , J )=0 . 

ERXCI, J)=0. 

ERY(I, J)=0. 

410 CONTINUE 
400 CONTINUE 

CxxxxxNUMERICAL INTEGRATION OF FEED POWERXxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DELTAP=PI/N5P 
DELTA=SI/N5P 
NSP1=NSPX . 5 
DO 201 1=1 , N5P1 
ZETA1=DELTA*(I-1) 

A145 = 2 . XQ 

FACT1=C05(ZETA1)XXA145 
FACT1=FACT1XSIN(ZETA1)XDELTA 
PINTN=PINTN+FACT1 
201 CONTINUE 

PINTN=PINTN/60 

PNUR=PINTN/PINT 

WRITE(9,993)PNUR 

993 FORMAT (/////, 5X, ’THIS NUMERICAL TO EXACT RATIO =*,F10.3) 


CxxxxxxsET SURFACE COORDINATESxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DO 10 J=1 , NSP 
DO 11 1=1, NSP 
X$=XSMAX-S5P*(I-1) 

YS=YSMAX-SSP*( J-l) 

RS=SQRT( (XS-XC)XX2+(YS-YC)XX2) 

IFCRS . LE. RMAX)GO TO 12 
GO TO 11 
12 CONTINUE 

Z$=-(XSXX2+YS#x2)/(4.*F)+F 

CxxxxxxCALCULATE THE NORMAL TO SURFACEXxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

AMAGN=SQRT(XSXX2+YSxx2+(4 .XFXX2)) 

ANX=-X$/AMAGN 

ANY=-YS/AMAGN 

ANZ=-2.XF/AMAGN 

CxxxxxxCALCULATE INCIDENT RAYxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

AIMAGN=SQRT( CXS-XF)xx2+(YS-YF)xx2+(ZS-ZF)xx2) 

SIX=(XS-XF)/AIMAGN 

SIY=(YS-YF)/AIMAGN 

SIZ=(ZS-ZF)/AIMAGN 

CxxxxxxCALCULATE REFLECTED RAYxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SRX=SIX-2.X(SIXXANX+SIYXANY+SIZXANZ)XANX 
SRY=SIY-2 . x C SIXx ANX+SIYXANY+SIZXANZ) XANY 
SRZ=SIZ-2.X(SIXXANX+SIYXANY+SIZXANZ)XANZ 
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0005700 

0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006700 

0006800 

0007000 

0007100 

0007200 

0007210 

0007300 

0007400 

0007500 

0007600 

0007700 

0007800 

0007810 

0007900 

0008000 

0008100 

0008200 

0008300 

0008500 

0008600 

0008612 

0008700 

0008710 

0008712 

0008800 

0008810 

0008820 

0008900 

0009000 

0009100 

0009200 

0009300 

0009400 

0009500 

0009600 

0009700 

0009800 

0009900 

0010000 

0010100 

0010200 

0010300 

0010400 

0010409 

0010410 

0010500 

0010600 

0010700 

0010800 

0010900 

0011000 

0011100 

0011200 

0011300 

0011301 

0011310 

0011315 


CxxxxxxAPERTURE COORDINATE OF REFLECTED RAYxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XA=X5“Z5* C SRX/SRZ) 

YA=YS“Z5* ( SRY/SRZ) 

CxxxxxxLOCATE THIS RAY IN THE APERTURE GRID MATRIXXXXXXXXXXXXXXXXXXXXX 


DO 21 N=1,NAP 
XCN)=XAMAX-SAPX(N-1) 

XDELTA=AB5CXCN)-XA) 

IFCXDELTA . LE. ALIM) GO TO 30 
21 CONTINUE 

GO TO 500 
30 CONTINUE 

DO 20 M = 1 t NAP 
Y(M)=YAMAX-SAPX(M-1) 

YDELTA=ABS(Y(M)-YA) 

IFCYDELTA.LE.ALIM) GO TO 40 
20 CONTINUE 

GO TO 500 
40 CONTINUE 

Cxxxxxx$TAR THE CALCULATION OF FIELDS AT APERTUREXxxxxxxxxxxxx 

A11=SQRT( (XC-XF)xx2+(YC-YF)*x2+(ZC-ZF)xx2) 

A12=SQRT( (XS-XF)XX2+(YS-YF)XX2+(ZS-ZF)XX2) 

ADOT=(XC-XF)x(XS-XF)+(YC-YF)X(YS-YF)+(ZC-ZF)x(ZS-ZF) 

AD0T=AD0T/(A11XA12) 

FPAT = ADQT*xq 

CxxxxjHIS IS TO CALCULATE THE TOTAL POWER INTERCEPTED BY REFL ECTORxxxx 
IFCFPAT.LT. IE-5) FPAT=lE-5 

CxxxxxFINISH FEED PATTERN CALCULATIONSxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
CxxxxxINCIDENT UNIT VECTOR CALCUL ATIONxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
EIX=- ( SIXXSIY) 

EIY=(SIX#x2+SIZ*x2) 

EIZ=~ ( SIZXSIY) 

A13=SQRT(EIXXEIX+EIYXEIY+EIZXEIZ) 

EIX=EIX/A13 
EIY=EI Y/Al 3 
EIZ=EIZ/A13 

Cx xxx x INCIDENT FIELD AMPLITUDExxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
EINX=EIXXFPAT/A12 
EINY=EIYXFPAT/A12 
EINZ=EIZXFPAT/A12 

CxxxxxREFLECTED FIELD AMPLITUDE CALCULATIONxxxxxxxxxxxxxxxxxxxxxxxxxx 
ERFX=2.X(ANXXEINX+ANYXEINY+ANZXEINZ)XANX-EINX 
ERFY=2.X(ANXXEINX+ANYXEINY+ANZXEINZ)XANY-EINY 
ERFZ = 2 . x ( ANXXEINX+ANYXEINY+ANZXEINZ) XANZ-EINZ 
CxxxxxxTHE PHASE CALCULATION xxxx* xxxxxxxxxxxxxxxx* *xxxxxxxx#*kx*x*kk** 
A31=SQRT( CXA-'XS)X(XA-XS) + (YA-YS)X(YA-YS)+ZSXZS) 
PHASE=C2XPI/AWAVE)X(A12+A31) 

ERX(N,M)=ERX(N,M)+ERFX 
ERY(N,M)=ERY(N,M)+ERFY 
PH(N,M)=PHCN,M)+PHASE 
PNUMCN,M)=PNUM(N,M)+1 
GO TO 11 

500 WRITEC 9^300) 

300 F0RMATC5X, 'CAN NOT LOCATE THIS RAY???????????????') 

11 CONTINUE 

10 CONTINUE 

CxxxxxTHIS WILL TAKE THE AVERAGE AMPLITUDE/PHASEXXXXXXXXXXXXXXXXXX 
DO 60 1=1, NAP 
DO 61 J = 1 , NAP 
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0011318 

0011320 

0011321 

0011322 

0011325 

0011326 

oomoo 

0011410 

0011417 

0011425 

0011431 

0011437 

0011500 

0011600 

0011700 

0011800 

0011820 

0011825 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012401 

0012500 

0012600 

0012700 

0012800 

0012900 

0013000 

0013100 

0013200 

0013300 

0013400 

0013500 

0013600 

0013700 

0013800 

0013900 

0014000 

0014100 

0014200 

0014300 

0014400 

0014500 

0014600 

0014700 

0014800 

0014900 

0015000 

0015100 

0015200 

0015300 

0015400 

0015500 

0015510 

0015511 

0015600 

0015700 

0015800 

0015900 

0016000 

0016100 

0016200 

0016300 

0016400 


IFCPNUMCI,J).EQ.O.) GO TO 61 
ERXCI,J)=ERXCI,J)/PNUMCI, J) 

ERYCI, J)=ERYCI, J)/PNUMCI, J) 

PHCI,J)=PHCI,J)/PNUMCI,J) 

61 CONTINUE 
60 CONTINUE 

CxxxxxxcOLLAPSE THE FIELD BY ROWS OR IN X, xxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITEC55)NAP,NAP 

WRITEC55)CCERYCI,J),J=1,NAP), 1=1, NAP) 
WRITEC55)CCPHCI,J),J=1,NAP),I=1,NAP) 

WRITEC55HC ERX( I»J),J=1,NAP) » 1=1 , NAP ) 

WRITEC55MCPHCI, J ) , J=1 ,NAP) , 1=1 ,NAP) 

CxxxxxxTO GET A COLUMN VECTOR OR Y DISTRIBUTIONxxxxxxxxxxxxxxxxxxxxxx 

DO 15 1=1, NAP 
DO 16 J=1 , NAP 


EXR(I) = EXRT(I) + ERXCJ,I)XCOS(PH(J,I)) 
EXI(I)=EXI(I)+ERX(J,I)XSIN(PH(J,I)) 
EYRCI>=EYRCI)+ERYCJ,I)xC0SCPHCJ,I)) 

EYICI)=EYICI)+ERYC J,I)XSINCPHCJ,I)) 

16 CONTINUE 

15 CONTINUE 

CxxxxxTHIS IS TOTAL POWER IN X AND Y COPONETS OF N . F. xxxxxxxxxxxxx 
CxxxxxxAPLYING THE FFT OR APERTURE INTEGRATION METHODXxxxxxxxxxxxxxxxxx 
NAP=NAPX5 

CALL FFT ( EXR >EXI#NAP»NAP»NAP,1) 

CALL FFT ( EYR, EYI , NAP, NAP, NAP , 1 ) 

A133=NAPXSAP/AWAVE 
NMAX=A133XSIN(ZETAXPI/180. ) 

NMAX=NMAX+1 
DO 70 1=1 , NMAX 
J1=NMAX- ( I — 1 ) 

A134=J1/A133 

ANGCI)=-ARSIN(A134)X180./PI 

J2=I+NMAX 

A135=(I-1)/A133 

ANGCJ2)=ARSIN(A135)X180/PI 

70 CONTINUE 

DO 71 J = 1 , NMAX 
I1=NAP-NMAX-1+1+J 

FIELDX(J)=FACTx(EXR(Il)XEXR(Il)+EXiai)XEXI(Il))xCOS(ANGCJ) 
XPI/180 )XX2 

FIELDYC J)=FACTX(EYR(I1)XEYR(I1)+EYI(I1)XEYI(I1) ) 
FIELDCJ)=FIELDX{J)+FIELDYCJ) 

12= J+NMAX 

FIELDX(I2)=FACTX(EXR( J)XEXR( J)+EXI(J)xEXI( J) )XC0S(ANG( J) 

XPI/1801XX2 

FIELDY(I2)=FACTX(EYR( J1XEYRC J)+EYI(J)XEYI( J) ) 
FIELDCI2)=FIELDXCI2)+FIELDY(I2) 

71 CONTINUE 
N1=2XNMAX 

CALL SCLBAKC .FALSE. , N1 , FIELD, RTNARR) 

GAIN=10 . XALOG10C4. XPIXRTNARR(2)XSAPXX4/PRAD) 

WRITEC 9 , 950 )GAIN 

950 FORMATC///// , 5X, 'THIS IS THE GAIN VALUE =’,E15.5) 

WRITEC 9 , 965)N1 

965 F0RMATC/////,5X, ’TOTAL NUMBER OF CALCULATED FAR-FIELD 

VALUES =»,I5) 

DO 73 1=1, N1 

FIELD(I)=FIELD(I)/ RTNARR (2) 

IF(FIELD(I).LT. 1E-6)FIELD(I)=1E _ 6 
FIELDCI)=10.xALOG10CFIELDCI)1 
WRITEC7,565)ANG(I) .FIELDCI) 

565 F0RMATC5X,E15.5,5X,E15.5) 

73 CONTINUE 

STOP 
END 
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Figure L- Reflector configuration. Frequency. 11.75 GHz. Polarization 
of feed: a-1 b-l/^Xand ^ • 90° for right-hand circular 
polarization; U E (0) - (cos 6) 3 - 6 ; U H (0) ■ (cos 0)^®. 
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Figure 2. - Reflector Z^ being illuminated by incident field from point source at F^. 
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